THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE. 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


*w_  i ^-— 


wmmm 


mad  \pruir«a&Tophnir*a|  infnrmafinn  Urranou 

llluU  OuillMCa  lUOIISHUUI  illlUIIIIUUUli  ilgUISUjf 

Because  of  our  limited  supply,  you  are  requested  to  return  this  copy  Vt'HEN  IT  HAS  SERVED 
YOUR  PURPOSE  so  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 
will  be  appreciated. 


ranCE:  when  government  or  other  DRAwIrfaS.  SPECIFICATIONS  OH  OX  HER  DATA 
iKEUSED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WTEH  A DEFINITELY  RELATED  | < 
JOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  S,  GQyERNMENT  THEREBY  INCUHo  \ \ 

fO  RESPONSIBIUTY,  NOR  ANY  OBLIGATION  ^'HATSOEVER;  AND  THE  FACT  THAT  THE  i 

GOVERNMENT  MAY  HAVE  FORMULATED.  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE  j 

SAID  DRAWINGS,  SPECIFICATIONS,  OR  CXrHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
EMPIJOATION  OR  OTHERWISE  AS  IN  ANY  MANNER  UCENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  FIGHTS  OR  PERMISSION  TO  MANUFACTURE,  • 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 

Reproduced  by 

DOCUMENT  SERVICE  CENTER 


KNOTT  BUILDING.  DiLYTON,  2,  OHIO 


k.. 

r ■ 

COLUMBIA  UNIVERSITY 
Hudson  Laboratories 
Dobbs  Ferry,  N.  Y. 


CU-31-54-0NR-271-Phys 


PROJECT  MICHAEL 
Contract  N6-0NR-27135 


W.  A.  Nierenberg 
Director 


Technical  Report  No,  24 


Analysis  of  a General  System 
for  the  Det<'c*.ion  of 
Anipli+'ii'^e  Lv.--.,uj.ated  Noise 


E*  Parzeii  aiio 


S.  Shiren 


Research  Sponsored  by 
Office  of  Naval  Research 


i 

S 

j 

i 

1 

I 

f 

I 

I 

i 

j 

j 

I 


August  2,  1954 


2 


TABLE  OF  CONTQn?S 


1. 

o 


3. 

4. 

5. 

6 1, 


n 

I ■ 


H , 

9 = 

10. 


Abstract 

Introduction 

Pov/er  spectrum  after  a square-lavi  detector 
The  detection  criterion 

Limiting  forms  for  the  detection  criterion 
Cross -correlation  before  the  multiplier 
Autocorrelation  after  the  multiplier 
Computation  of  the  detection  criterion 
Application  to  the  case  of  Gaussian  filters 
Extension  to  the  inclusion  of  background  noise 
Conclusions 

Appendix:  The  mathematics  of  noise 

References 


LIST  OF  FIGURES 


Figure  1 
P'igure  2 


General  detection  system 
A typical  detection  scheme 
modulated  noise 


‘or  amplitude 


- '-to 


3 

4 
A 


9 

12 

15 

21 

23 


26 

oo 


33 

d4 


37 


Page 

38 

39 


“r*'- 


m 


I 


!? 

i 

Jt 

.55 


. i 

15 


I 


BSTRACT 


A general  system  for  the  detection  of  amplitude-modulated 
noirie,  in  the  presence  of  background  noise,  is  analyzed  with  a 
view  tovjard  determining  the  behavior  and  optimum  design  of  the 
system.  The  unmodulated  noise  carrier  and  the  background  noise 
arc  assumed  to  be  independent  eta  ionary  Gaussian  random  time 
fur.;- cions.  The  modulating  functlo:  is  a time  function  which, 
in  general,  is  random,  non-Gaussian , and  nonstationary.  A 
detection  criterion,  as  a moasur..  the  performance  of  the  sys- 
\>em,  i.'  defined,  and  computed  in  terms  of  the  input  power  spectra 
and  the  transfer  characteristics  of  the  system.  The  techniques 
used,  and  the  intermediate  mathematical  results  obtained,  ar*^  of 
interest  in  themselves.  The  results  are  applied  to  give  a 
detailed  analysis  of  a typical  d3toction  system  which  is  a 
special  case  of  the  general  sy  t_Ji. 
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1.  INTRODITCTION 

The  problem  of  the  detection  of  signals  having  the  character 
of  amplitude  modulated  random  noise  is  of  recent  interest  in 
communications  theory.  For  example j a paper  by  Deutsch^  treats 
the  effect  of  a linear  detector  on  such  signals  when  the  modula- 
tion is  periodic. 

In  this  report,  we  treat  a general  detection  system  involv- 
ing square-law  detectors.  We  shall  compute  various  statistics  of 
the  output,  with  and  without  modulation,  and,  from  these,  compute 
a detection  criterion  which  we  define  as  a measure  of  the  detecta- 
bility of  the  modulated  noise. 

The  system  is  diagramed  in  Pig.  1.  The  input  filters  and 
include  the  characteristics  of  any  receiving  equipment  which 
might  ordinarily  precede  them.  The  only  assumptions  made  on  the 
filters  and  are  that  they  do  not  pass  dc.  The  output  inte- 
grator is  perfect  with  integration  time  T. 

We  take  the  input  u(t)  to  be  a stationary  Gaussian  noise  y(t) 
modulated  with  index  m by  a signal  inunction  g(t),  statistically 
independent  of  y(t).  Thus 

u(t)  ry(t)  j^l4-mg(t)J  . (1.1) 

We  assume  both  the  noise  and  the  signal  to  have  zero  mean  level; 
that  is*, 


<yCt)>  = 0,  s:0  . 

(1.2) 

Then,  the  total  input  power,  P^,  is 

^u  ~ ^)t  Av  ”* **n 

(1-l-m^Pg) 

(1.3) 

where 

Pn=<y^(t)>  =1% 

( CO ) dci> 

(1.4) 

; (cO  ) du) 
g 

0 


* This  notation  is  defined  in  the  Appendix  in  which  we  summarize 
the  main  notations  and  notions  on  random  time  functions  that  we 


use. 


Thus  m~?  is  the  ratio  of  the  sideband  power  to  the  carrier  power. 
S 


We  shall  assume 


r>  c'' 


since  this  is  the  interesting  case  in  a detection  problem.  The 
T?ethods  usedj  however j are  applicable  for  any  value  of  the  ratio.- 

If  by  detection  is  meant  the  extraction  of  information  about 
the  modulating  function  (which  would  also  seem  to  be  the  best 
means  for  distinguishing  this  type  of  signal  from  unmodulated 
noise),  then,  the  inherent  presence  of  fluctuation  noise  at  the 
output  of  the  detector,  due  to  demodulation  of  the  input  noise 
carrier,  makes  this  a problem  of  the  detection  of  signals  in  the 
presence  of  additive  noise. 

In  order  to  see  this  more  clearly,  and  to  better  understand 
the  action  of  the  general  system  of  Fig.  1,  we  shall  compute  the 
output  power  spectrum  of  a square-law  detector  when  the  input  is 
unmodulated  noise.  This  calculation  also  gives  a simple  illustra- 
tion of  the  methods  used  in  treating  the  general  system. 
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2.  PO\m  SPECTRUM  AFTER  A SQUARE-LAW  DETECTOR  | 

1 

? 

If  the  input  of  a square-law  detector  is  given  by  (1.1), 
then,  the  output  /(t)  of  the  squarer  is  given  by  i 

v(t)=  u^(t)  - y^(t)  [l+2mg(t)  -h  mV(t)  ] . (2.1)  | 

In  order  to  cornpute  the  power  spectruim  G (n)  ) of  v(t),  we  | 

first  compute  the  autocorrelation  R^(t  ),  and  then  obtain  ) | 

by  means  of  the  Wiener-Khintchine  theorem  (Eq.  (Ap)  in  the  Appendix).  I 

R^,(r  ) :=  ^v(t)  v(t  -+T  ) ^ (2.2)  1 

= {y^(t)  y^(t  + 't  )')  I 

+ (g^(t + t ))  4-  4 <Jg(t)  g(-b+^))  i 

4-  m^  /g^(t)  g^(t  t-n; ) ^)  I . I 

The  terms  in  (2.2)  which  involve  m and  m^  vanish,  since  we  ] 

assume  that  all  time-ensemble  averages  of  odd  powers  of  g(t)  vanish.  I 

To  compute  the  foregoing  fourth  order  ensemble  average  of  y(t).,  | 

we  use  the  formula  (A7)  given  in  the  Appendix,  which  expresses  the  I 

higher  order  moments  of  a Gaussian  random  time  function  in  terms  of  s 

its  autocorrelation  function  R^v  V ) . | 

^ I 

We  have  • 

-'2.  ! 

(y^(t)  y^(t+t))=(y^(t))  {y^(t -ht  ))-f  2 [{y(t)  y(t  H-T:' ))]  (2.3)  \ 

2 " I 

- Ry^(O)  4-  2 ).  j 


As  regards  the  fourth  order  tine-ensemble  average  of  g(t),  we 
make  the  assumption  (which  is  fulfilled  by  a large  numoer  of  random 
time  functions)  that 


(s^it)  j^^=Rg^  (0)4-f(T) 


(2.4) 


where  f(X  ) is  of  the  order  no  greater  than  that  of  R (0). 

o 
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For  example,  if  g(t)  has  Gaussian  statistics,  then 

f(r  ) = 2 (r)  ^ 2 Rg^  (0). 

Consequently,  we  obtain  from  (2.2)  that 

R^(Tf  ) =r  [p.y-  (0)  + 2 Ry^  (-t  )}  ([l+  R^  (0)  J ^ 

+ 4 R„  (r  ) S-  m'’’  )} 

B J * 


(2.5) 


(2.6) 


In  view  of  ass^jmption  (1.5)»  and  the  foregoing  remarks  about- 
f('t'),  we  may  henceforth  ignore  the  term  in  (2.6)  involving  m^. 

The  power  spectrimi  G^(0)  is  now  easily  obtained  by  taking 
the  Fourier  integral  of  (2.6).  The  resulting  terms  group  themselves 
naturally  into  three  groups;  by  (1.4),  we  write  P and  P for  R-(0) 

IX  g y 

and  R^(0)  respectively. 

o 


The  power  spectrum  G (w)  is  the  sum  of  dc  terms:  (2.7) 

[l  + m^  Pg]  ^ , 

2 2 

steady  state  or  signal  terms:  4m  P^  Gg(W  ) , 

noise  or  fluctuation  terms*  ^1+ m^  ^g  ^ ^ Gy(u>->j.  ) 

2 Jd-v  G (v  ) J d^  G (^jj.+  w-V  ). 


The  integrated  power  is  then,  dropping  terms  involving  m Pg  by 
assumption  (1.5)> 


dc 

P ^ 
n 

signal 

4 m' 

noise 

2 P 

- , 2 


■n  g 


n 


(2.8) 


Thus,  as  stated  at  the  end  of  Sec.  1,  at  the  output  of  the 

square-law  detector  the  problem  reduces  to  that  of  detecting  a 

2 

signal  of  total  power  4 m presence  of  an  additive 

noise  of  total  power  ‘xi* 

This  problem  is  usually  attacked  by  analyzing  the  detector 
output  with  a band-pass  or  low-pass  filter  followed  by  squaring 
and  averaging  to  determine  the  mean  square  signal  and  noise. 

It  is  easily  seen  that  the  complete  detection  system,  ; 

diagramed  in  Fig.  2,  consisting  of  receiver,  square-law  detector,  ] 

-.1 

filter,  squarer,  and  finite  time  integrator  is  a special  case  of  | 

the  system  of  Pig.  1,  obtained  by  setting  H » = Hg  and  K »=•  * Kg*  j 

This  case  is  treated  in  detail  in  Sec.  8.  I 


I 


I 
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3.  THE  DIRECTION  CRITERION 

Suppose  that  it  is  possible  to  measure  the  output  w(T)  of 
the  system  of  Fig.  1 when  the  input  is  unmodulated  noise.  A large 
niamber  of  such  measurements  will  have  an  average  and  a mean  square 
fluctuation  given  by  the  ensemble  average  ^w(T)^  ^ and  the  variance 

A 2 

• t W~fT 

n 


'J 


(3-1) 


Similarly,  a large  number  of  measurements  of  the  output  when 
the  input  is  amplitude  modulated  noise  will!  have  mean  value  ( w(T)), 
and  variance 

2 _ /-2  2 


mn 


0" 


mn 


- - <*«)> 


mn 


(3.2) 


The  subscript  n on  a quantity  is  to  indicate  that  the  quantity 
is  to  be  evaluated  under  the  assimption  that  the  input  to  the  system 
is  unmodulated  noise,  while  the  subscript  mn  will  be  used  to  indi- 
cate that  the  quantity  is  to  be  evaulated  under  the  assumption  that 
the  input  is  amplitude  modulated  noise. 

Under  the  assiunption  that  the  observations  are  independent, 
the  difference  A(T)  of  the  observed  outputs  in  the  two  cases, 


A(T)  = w(T) 


mn 


w(T), 


(3.3) 


will  have  mean  value  \A(T)^  =■  ^ w(T)^  and 


V MX  X.  <W  ^ M.  V ^ ^ 

mn  u 


Using  ■th.^ory 

ing  hypot.V.eses , if  the  probability  distribution  of  A(T)  is  known, 
by  testing  the  hypothesis  that  ^A(T)^  = 0 one  can  detect  the 
presence  or  the  absence  of  signal.  To  any  preassigned  probability  p 
one  can  find  a number,  K,  such  that 

iA(i)i  i K (3.4) 

with  probability  p.  As  the  detection  process,  one  then  adopts  the 


1C  - 


following  rule.  One  decides  that  noise  alone  is  present  if  the 
difference  of  the  observed  outputs  lies  in  the  range  given  by 
Eq.  (3»4-),  and  that  signal  (in  our  case,  modulated  noise)  is 
present  if  the  difference  of  the  observed  outputs  lies  outside 
the  range  given  oy  Eq,  (3*4-). 

Motivated  by  the  foregoing  considerations,  we  therefore 
adopt  as  o\xr  detection  criterion,  denoted  by  D(T),  the  following 


ratio  * 


D(T) 


L wir. 


^ 1 1 
‘^n  j 


(3.5) 


It  is  presumed  that  the  larger  D(T)  is  for  a given  detection  sys- 
tem, the  better  will  be  the  system's  performance,  and  the  more 
likely  it  will  be  to  detect  the  presence  of  signal. 

In  many  cases,  it  can  be  assumed  that  the  fluctuation  term 
for  the  case  of  signal  and  noise  in  the  input  is  roughly  equal  to 

the  fluctuation  term  01  for  the  case  of  noise  alone  in  the  input. 

n 

We  then  have  for  the  detection  criterion 

])  (T)  = 


% — >•  • / / 


JJL. 


n 


(3.6) 


In  view  of  assumption  (1,5)  j one  could  show  that  is 

roughly  equal  to  <r_,,  and  we  therefore  use  Eq.  (3.6)  for  the 
detection  criterion. 

We  now  turn  our  attention  to  the  problem  of  expressing  D(T)  in 
terms  of  the  statistics  of  the  noise  and  of  the  modulating  function. 
This  problem  is  solved  in  several  stages,  as  follows. 

In  Sec,  4,  we  express  D(T)  for  large  T in  terms  of  the  cross- 
correlation, ^ ( TT ) , of  the  outputs,  v-^(t)  and  V2(t),  after  the 
second  filters,  and  the  autocorrelation,  R^('T?),  of  the  output,  w(t), 
of  the  multiplier.  We  also  obtain  a similar,  but  more  complicated, 
expression  for  D(T)  for  small  T. 
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In  Sec.  5>  the  cross-correlation,  p ),  is  computed  under 
two  different  assumptions.  First,  it  is  computed  under  the 
assTomption  that  the  input  of  the  system  is  unmodulated  stationary 
Gaussian  noise,  in  which  case  we  denote  it  by  Since 

7 •'In 

bhe  stationary  character  of  Gaussian  noise  is  preserved  after 
passage  through  non-linear  devices,  it  holds  that 

+ Vtav]  n 0 = 7) 


= ^Vi(t)  Vg(t  -f- 


Second,  p('T^)  is  computed  under  the  assumption  that  the  input, 
u(t),  is  amplitude-modulated  noise.  In  this  case,  we  denote  it 

Thus , 


"^2^*  >>TAv] 


mn 


(3.8) 


In  Sec.  6,  the  autocorrelation,  Is  computed  under 

the  assumption  that  the  input  is  noise.  We  have  therefore 


^ = \ <(*00  +■«  )>TAv  ] n 


(3.9) 


— (^w(t)  w(t  X ) \ 

% / ** 


Finally,  in  Sec.  y,  the  results  of  the  preceding  sections  are 
combined  to  yield  the  required  expression  for  D(T). 
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4.  LIMITING  FORMS  FOP.  THE  DETECTION  CRITERION 


A Formula  for  D*^(T) 


Since 

fT 

^?(T)  =1  w(t)dt 
Jo 


it  folloTis  that 


(t)  V2(t)dt, 


J O 


mn  J 


{w^(T))  f r dt  dt'  /w(t)  ^(t')*)^ 

- o Jo 


(4.1) 


(4.2) 


(4.3) 


(4.4) 


r T r T fT 

= 1 dt  dt‘  R«(t'-t)  = 2j  drR^('V)  (T-Tf). 

Jo  Jo  >'0 


Therefore, 


2 PT 

(T))n  ' ” 2 at  R^j(r  ) (T-tr  ) (4.5) 

2 

- [Tfp  ^0)]  = > -Pn  <°0 


P 


Corsequently  we  may  write  for  D Vi’Js 

r . rT  . 


?^nr)  ^ 

T z fO-^/r)  [P-M- io)] 


(4.6) 


I 
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The  Limit  for  T Large 


If  thfi  noise,  y(t)^  has  a continuous  power  spectrum 
as  we  assume  to  he  the  case,  then  it  will  be  shewn  in  Sec.  6 that 
the  power  spectrum  G^(co)  of  the  output  of  the  multiplier  (when 
the  input  is  noise)  will  also  be  continuous  except  for  a dc  delta- 
function  term  equal  to  2 (0).  Ytfe  therefore  define 

to  denote  the  continuous  power  spectrum  of  the  output  w(t)  of  the 
multiplier • 

The  autocorrelation,  corresponding  to  G’  ( <*>)  is  then 

given  by 


(0) 


(4.8) 


R*  (U  ),  unlike  R^('^),  has  a finite  integral  from  0 to^®  . By 
the  Wiener-Khintchine  theorem,  it  follows  that 

(r)  J-r  = 


Let  us  now  pass  to  the  limit,  as  T o®  > in  (4.6).  By  defini- 
tion, the  first  term  of  the  numerator  tends  to  (0).  From  the 

• I mn 

fact  that  R'^(Tf  ) is  integrable  it  follows  that  the  denominator 
tends  to  Tf  (0) . 

VI 

We  therefore  ootain  thai’  - 


I 


T TT 


/i  /rti  1 


7T  (o) 


(4ol0) 
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(4.11) 


Thus  for  T large,  we  have  approximately 

3)  (T^  = Jf  • (4.11) 

Thus  the  detection  ratio  increases  with  the  square  root  of  the 
integration  time.  In  other  words,  in  order  to  double  the  detecta- 
bility of  a weak  signs.l  in  the  presence  of  noise,  it  Is  neces.sary 
to  quadruple  the  integration  time. 


The  Limit  for  T Small 

For  T small,  the  denominator  of  Eq.  (4.6)  is  approxi.mately 
T “ p j * However,  there  is  no  simple  expression 

for  the  first  term  of  the  numerator,  unless  the  signal  is  station- 
ary, when  the  integration  is  no  longer  necessary.  We  therefore 
can  only  obtain  the  formal  expression 


2,  P*(T)  = f-i^o  T . 

Rw(0)  - p/(o) 


In  Sec.  5>  we  will  show  how  the  first  term  in  the  numerator  of  Eq, 

(4.12)  may  be  evaluated  if  we  possess  a knowledge  of  the  complete 
statistics  of  g(t);  that  is,  if  we  know  the  mean  value  function 


Ag(t)  = {^g(t)^ 


and  the  covariance  function 


r*  (tj^,  g(^2^^  * 


(4.13) 


(4.14) 


It  may  be  noted  that  the  assumptionr.  we  have  already  made  on  g(t)  may 
be  written 


0 = [ 

^ ^ J r Av 


(4.15) 


(4.16) 


, tMCTroiffit— 


sti^rnmuMi  iw-^u  ♦ly  # 


J ^lAl^mf:-^  's-'. 
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5.  CROSS-CORRELATIOK  BEFORE  THE  MULTIPLIER 


In  this  section,  we  compute  P (ly  ) and  P ) • We  make 
’ in  I mn 

the  convention  that  all  integrals  which  are  taken  from  - oo  to 
Of*  are  to  be  written  without  the  limits  of  integration.  Since 


a/i+r-A^7)u.(fc«-/V^y5* 


2) 


it  follows  that 


P 1 1'\  - 


\ X rr  I , I \ 


(5.3) 


= J'  <F;tstfc^3><F:[4-+r«3t^r)]>^^^  (5.4) 


where  for  brevity  we  use  the  single  primed  integral  sign  \ to  denote 
the  six  fold  integration, 


f 


and  for  a function  f(t),  we  define  the  notation  F,  [f(t)]  by 

FXW^- -f  (^  + 't'-fVy)  (5.6) 


In  Eq.  (5. 4),  we  have  used  the  formula 

< KtMrt'l  f.  i Ch ‘*>3>  <5-7) 

J -•  Triv  •'  t ny 


which  follows  frotD  the  statistical  independence  of  y(t)  and  g(t).  and 
the  stationary  nature  of  y(t). 


'-J-. 


4 

I 


\ 

3 


-31 

'4 


a 


,73 


J 

)J 


> . 


We  henceforth  drop  the  subscript  y in  writing  the  autocorrela- 
tion and  power  spectruu  of  y(t).  Since  y(t)  is  Gaussian,  we  may 
express  the  four  fold  ensemble  average  C y(t)3^in  terms  of 

the  autocorrelation  R("J*  ) of  y(t)  by  means  of  Eq.  (A7)  given  in 
the  Appendix.  We  thus  obtain 

{FltjU)])  "-=  (5.8) 

+ R ('t-h  ^ - '^1  R (t -Ho(2 


1-f  mg(t)J 

obtain 

F,  [l+wj 


is  easily  expanaea 


(5.9) 
r 

-f  I ^ ("t-  c^r  ^ ^ 5 2 t - “'‘5)  2 (i-i r 

^Tn^l  o<2.~'?)  ^(■b+x-(3,"nV 

+ 3 U+'?- (^rS)  T 

+ 3 U-  3 (-t+i'  3 (-t+x-  J 

T)n^j  q('b-^.-Vi  q s p,-‘n)  q 

d * Cl  C/  '• 
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Upon  taking  the  tlme-ensem'ble  average,  we  obtain 


<F;ti+i*i56tO>^^=  1 (5.10) 


{% 


ST* 


11 


The  terms  in  m and  m-'  vanish,  because  it  has  been  assumed  that 
all  time-ensemble  averages  of  odd  powers  of  g(t)  vanish.  In  what 
follows,  we  may  also  ignore  the  term  in  Eq.  (5»10)  involving  m^, 
since  we  assume  that  the  fourth  order  time-ensemble  average  of  g(t) 
is  of  the  order  of  (0).  Consequently,  the  term  in  Eq.  (5*8) 
involving  is  of  the  order  of  this  term  may  be 

dropped  in  view  of  the  discussion  in  Sec.  1. 

In  view  of  Eqs.  (5*10),  (5»8),  (5»4),  and  (5*3) » we  could  now 
write  an  expression  for  «hd  for  p^„( tf ) in  terms  of  the 

modulating  index,  m,  the  autocorrelation  functions,  R(t;)  and 
R (T),  and  the  filter  impulse  functions,  h^(ot),  hgC^),  )> 

k2<>|).  However,  it  is  more  convenient  to  express  the  cross- 
correlation,  p {'t  ) i in  terms  of  the  power  spectra,  G(c<))  and 
G (cd),  and  tne  filter  transfer  functions, 

and  To  do  this,  Eqs.  (5«10)  and  (5*8)  are  substituted 

in  Eqs.  (5.4)  and  (5«3)»  In  the  resulting  expression,  R('l' ) and 


R, 


g 


(t  ) are  replaced  by  the  Fourier  integrals  which  relate  them  to 


Qiu})  and  )»  given  by  2q.  (A6)  of  the  Appendix.  By  inter- 

changing these  Fourier  integrals  with  the  integrals  indicated  in 
Eq.  (5.5)?  and  performing  the  latter  integrations,  we  obtain  ths 
following  expressions  for  and  where  an  asterisk,*, 

denotes  a complex  conjugate.  In  deriving  these  formulas,  we  have 


I 
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made  use  of  the  facts  that  a spectral  density  Is  an  even  function 
of  its  argument;  while  any  filter  transfer  function,  say  H,  is 
Hermitian;  that  i^  H(-0)  ) c:  H ♦ ((*>)• 


= 

Ki(0)  K^to)  GiW 

+ 1 GW  ^(«»)  HtW 

K|  ('J.+W  /wj+'^iV 


(5.11) 


(5.12) 


K,{0)  Ki(o')  { I-U,  CiW  ^ 

~Tn"  1B,WP  jAvfijfv) 

ih’’  iclwi  G(W  IHjWI^Iav 

+ a Qk^  fik")  Hfw 


4 

+■ 


4 IV  kiWKaV; 

f Au,  gU)  II, (W  Ui('»ii+'*)  ^<*“1  1^2 W >') 


+ W-  IAv  S(-*)  flAwJ.  Awx 


H.(w.+v^  K',(w.+0  I^z’^(h+0 

+ flir  cJvL'i  sW  GW) 


H /u)^  W,*”  W-v)  U.(u")  H,*‘(w,+v)  fC.W+W  ki 
ai«"  la,  6j(v)  e'-'^  6(W 

■ i i I \ k ^ \ *_L  ^ k ■ • ' • - i - • \ \i  ^ i z .\  \.  \/  /*■»  1^.  ^ 1*1 . ii  jL 

f2V  I'Av  e"-'”  6(W 

H.(ui,i-t)  H,(W  K,W+*>.+v)  ifz  (u),+^i4->) 


TT 


•s.'jiVa.  **.**.,».  ^ 


VMnBMMHPPmiM 


hv 


R 


rW?H<HAfc.-  A-  ■A’, 
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The  reader  should  observe  that  all  the  foregoing  integrals 
are  real  valued  quantities. 

Later^  we  will  assunia  that  K-,  (0)  0 = KgCO).  Therefore,  we 

have  written  the  roregolng  expressions  for  pnm- ^ ' 

in  such  a way  as  to  make  evident  the  form  they  assume  in  this 
case. 


To  conclude  this  section,  we  indicate  how,  with  a knowledge 
of  the  mean  value  function  Eq.  (4.13)  and  the  covariance  function 
Eq.  (4.14)  of  g(t).;  one  is  able  to  evaluate  the  quantity 


T->0 


I 

'T 


< viW> 

Tnyi 


(?.13) 


required  in  Eq.  (4.12).  Using  the  methods  of  this  section,  it  is 
evident  that  Eq.  (5«13)  is  equal  to 

I <F;hw]>  JtotJ 


A 


Using  the  expansion  of  f 14-mg(t)] 
easily  seen  that,  up  to  terms  in  m^, 


given  by  Eq.  (?.9)>  it  is 


^ (■*•)]  > = i.  (5.15) 

4.  T-ft-’i)  + Pj  \ r - /*t-'n) 

+ Tq,  T- (i.-i)  4 r,  t:- i3,-n)  i 


53f 


KrTrTTjM. 

yfgmfe  wimietmwTiir^:  y>  r» 


6. 
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AUTOCORRELA.TXOK  AJFTER  THE  MULTIPLIER 


In  this  section,  v?e  compute 
/ ^ / 


vX 


Rur  VT‘>tt)  Vi(fcvxV- 


n. 


(6.1) 


the  autocorrelation  of  the  output  of  the  multiplier  vjhen  the  input 
is  unmodulated  stationary  Gaussian  noise.  Using  the  same  approach 
as  in  the  previous  section,  we  write 


RurW=  f 


(6.2) 


I ’ 

where  for  brevity  we  use  the  aouble-primed  integral  sign  J to  denote 
the  twelve  fold  integration 


jIII 

JjJJ  (^o(,  (iefj  Jo(j  ^r\|  (c<,j  (•^3)  ) 

fjji  4(W  W) 


(6.3) 


'ki 

I 


4 


and  we  define  the  notation  mean 

(i  - <J,  - '^,)  (-t-  o(j-X')  ij  (t  + 1-  /j-\)  I3  (■!:+•«-  »<,''?,)  (6.4) 

( t - (i, -•>),)  ^ C-t  - Pi-'li)  'j  ^ 


Now 


y(t)  ^^is  an  eight  fold  average  of  a Gaussian  random 


function.  Consequently,  by  Eq.  (A7)  of  the  Appendix,  it  may  be 
expressed  as  tiie  sum  of  105  terms,  each  term  a product  of  four 


4 

■s 

' a 


‘i 

I 

■ta 


a 

gs 

i 
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two  fold  averages.  It  is  easy,  although  tedious,  to  enumerate 
these  terms.  It  Is  found  that,  if  we  assume  that 


K. (0) = 0 = K2CO), 


then  4-5  of  these  terms  vanish  when  the  integrations  indicated  "by 
Eq.  (6.3)  are  performed.  By  the  same  method  as  used  in  the  pre- 
vious section,  the  remaining  60  terms  can  he  expressed  in  terms 
of  the  power  spectra  G(c»J  ) and  Gg(co),  and  the  filter  transfer 
functions  H-,  (co  ),  K^(a>),  and  K^Cto).  We  finally  obtain 

the  following  expression  for 

!v=  1?^.  (t)  = 

f 


i£, 

no 


^ 1 i i Hi"' I'a)  Hi(wj)  HzM 


rrrr 


If 

n 


+ M Jjjj  tW,  ^ 

H,w  Hj  (u.,') 

4-^t  IjI?  6(u>a)  6(1^3)  ^ I ' 

Hl(tw3Vr\z 

(uJi-Hi)a)  K jl  isi  Kz.  ( ^^au""  ^3^ 

-hi. ip  jj[]  ito,  JlvOi  iuJ)  ^(‘4  GU.)  ^ 

Hi(tjii)  Hu^wx)  Hi (^*^3)  H{(u)m)  Hn 

Ki  (lUi  + 5*)3')  K-1^  (u)| +•<«'»<)  Ki 

AC  A)  ^k) 

!Hz(w^Vl“  Hi(uj3)WzrK)  j-ir^W 
Ki(‘0jtW3l  K-^ 


\.^ 

-t 

-i 

iT 
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■■.« 
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7.  CCaiPUTATION  OF  THE  DETECTION  CRITERION 


In  vievi  of  Eq.  (4.11),  the  detection  criterion  D(T)  is 
known  for  T large  as  soon  as  we  know 


I-"  , (0) 

[mu' 


.O  ' — 


(7.1) 


2-JT  G'^(O) 


(7.2) 


These  quantities  inay  now  be  expressed  in  terms  of  the  charac- 
teristics of  the  filters  and  the  statistics  of  the  input  noise  and 
modulating  function. 

From  Eqs.  (5.11)  and  (5.12),  and  the  assumption  Eq.  (6.5),  it 
follows  that 


(7.3) 

^7:{^  ]jL>6^(v)  K^(vj  (q(w,) 

GiU')  Hi"" 

jdl'ii  ^^(v)  fo+^)  (f^i) 

Hi  Kj  (W|4-u>0) 

V Jctv  jj  ^U)  HiW 

r4i(*<J2^  Hx  (wa.*h>^  iv>.^ 

i'iv  o^(y)  fj  <iu)|  Hi 

U ( n I ^,  ( U).  ^ w*.  ^y)  Kx*  (u),  + u>x + v) 

4 i Trr  TjLv  Hi(u),^-v] 
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It  is  immediately  seen  from  Eq.  (5«H)  that,  under  the  assump- 
tion Eq.  (6o5),  the  first  term  in  the  expression  Eq.  (6.6)  for 
is  equal  to  Therefore  is  given  hy  the 

r^her  terras  of  Eq.  (6.6).  By  integrating  these  terms  -with 
respect  to  't  from  - to  obtain  the  following  express- 


Zn  Gur(O)  = ^TT 

Uick1>)  ffd^  hMhJ^o^ 

H2.  J 

4-  j’^/i  I a(ju,)l^  //<i(h)/^  IHz(m>I^  ^ 

Ifdu  Hz^(u)tA)! 

4-  }d^  I iwl  ,2 

{/fjlw  A',(w)  AiVw)  6f<^+vw) 

fdM  ffdu;,  Qk)€(i^^V^) 

/a 

Aj  lc0,-hc>^2 K2* 


BSSS 


- 25  - 


From  the  above  results,  ve  may  immediately  obtain  the 

detection  criterion  D(T)  for  the  system  of  Fig.  2 by  setting 

Ht(W  ) HgCtO  ) ^H(w)  and  ^ ).  D(T)  for 

J.  ^ . /A  TT  % vi<-\mu  -hVio  +.PTms  on  the 


i.11  \WV  i ..rj  - a.  1-U_ 

^ , vtr  (A.  IT  ^ hut  no™  t.hfi  terms  on  the 

T large  is  stxiJ-  givau  by  ^.q.  , euu  u«-..  — - 

right  hand  side  of  Hq . (4-, 11)  may  be  expressed  as  folloiNSJ 


(7.5) 


jkv  G h)  IK(y)r  lU  Giu)  H*(j) 


-J-Xth’-  fiy  Gi.Ci)  Iliwi  G(‘^^  //'/ftv) 


..H  u 


/ (4 1 i),  •t-y')  I ^ 11^  i 


4-y  *,»■  fjv  ^(v)  ffju,  GM  6fe) 


2,rCA')= 

ix  Ua  /'<(>)!''  ( u-  ^ 

+ 2 /ci^  166*)!^  kut  G(e^A^  } 

-j-  IdixL  (k'W!^  fUi>>a^z  g0->)  gM  ) 

/rt^K)P  )H^(4+a)/^  J 


1 

4-  <,4. 


& 


- 26  - 


8.  APPLICATION  TO  THE  CASE  OF  GAUSSIAN  FILTERS 


To  illu.strate  these  results ^ let  us  compute  D(T)  for  the 
system  of  Fig.  2.  For  mathematical  convenience,  let  us  assume 
that,  up  to  phase  factors,  the  filter  transfer  functions  H(cJ  ) 
and  K(0  ) are  given  by  Gaussian  functions,  as  follows: 


, , r ; /U-J1..V1  [-  I / ui  1-A^Y  1 


IC  M = 1 5 * m.  t T (^^Tl 


(8.1) 


(6.2) 


While  not  physically  realizable,  Gaussian  snapped  filters  are 
often  good  approximations,  for  mathematical  purposes,  of  actual 
filters.  Under  this  assumption,  the  evaluation  of  the  integrals 
in  Eqs.  (7*5)  and  (7*6)  is  much  simpler  than  it  would  be  otherwise, 
for  we  may  use  the  following  useful  formula  for  the  product  of  two 
Gaussian  factors: 


jW 


r > ( AV  1 


= ‘I'fft-i  J -M<f  L-t 


err-  t- 


(8.3) 


where 


il  = 


07  

>**•  ^ i_  — ^ 
^ ^ X 


(8.4) 


« -i : 
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It  is  ea.^y  to  verify  Eq.  (8.3)  by  expanding  both  sides, 
Eq.  (8.3)  it  folloTus  that 

Jj. 


From 


(8.5) 


For  the  signal  pov.er  spectrum  G^(<«>),  we  will  assume  a 
Gaussian  function  centered  at  roughly  the  same  point  as  is  K(ui). 


r -L  cr;'^  o/to  [-■!- 


' I ^ i \ Oj  i J 

f l / tb+  1 T 

4 ^L-1  J J 


(8,6) 


There  is  no  difficulty,  however,  in  treating  any  other  form  of  sig- 
nal spectrum. 

The  noise  power  spectrum  G(cb;i  is  assumed  to  be  flat  and  iden- 
tically equal  to  a constant  G. 

We  assvune  that  all  crOss-product  terms  may  be  ignored  when 
these  Gaussian  functions  are  multiplied.  We  then  have  by  Eq.  (8.3) 


k V 

MM  M(“+M  = 
r r 1 [^-AjrM.hY  1 , r. 


(8.7) 

_ i 1 1 

. ^ tJ/2  JJ 
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We  also  sssume  that 
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k 

X 


(6.8) 


1^ 


There  Is  now  no  difficulty  In  evaluating  Eqs.  (7«5)  and  (7«8). 
We  obtain  the  following  approximate  expressions: 


pm.  I o')  - P.^0)  = 

* r.  J • ii  J 


HtT  7« 


2ir  G'ur(O)  = 


9*^  I /-TW  I 


^o 


i _ 


^it 


r:- 


{i+ 


(8.9) 


(8.10) 


i 


u 

Vi 


ks  a measure  of  the  ratio  of  the  bandwidth  of  the  spectrum  of 
the  modulating  function  to  the  bandwidth  of  the  power  transfer 
functions  and  lK(a>)| 

respectively, 


of  the  filters  H and  K,  define, 


R. 


:=  ^ 
(T^ 


il  . 


(3.11) 
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9.  EXTENSION  TO  THE  INCLUSION  OF  BACKGROUND  NOISE 


It  is  of  sorr.a  interest  to  consider  the  problem  of  detecting 

the  amplitude  modulated  noise  in  the  presence  of  background  noise. - 

We  let  z(t)  denote  a stationary  Gaussian  random  time  function, 

Yiith  autocorrelation  R_(T")  and  power  spectrum  G^(Ci>).  Vi/e  assume 

z z 

z(t)  to  be  statistically  independent  of  y(t)  and  g(t).  Let 


(t)  iT.  y(t)  4*  z(t) . 


(9.1) 


Then  y^^Ct)  is  a stationary  Gaussian  random  time  function  with 
autocorrelation 


R]  ( ^' ) — ) -t"  ^ ' 


/a  2' 


and  power  spectrum 

G^(o>)  n Gy(CJ)4- 


(9.3) 


If  the  input  of  the  system  of  Fig.  1 is  y^(t),  then  the  cross- 
correlation  of  the  outputs  before  the  multiplier  (denoted  by 
and  the  autocorrelation  of  the  output  after  the  multiplier,  continue 
to  be  given  by  Eqs.  (5.11)  and  (6.6),  respectively,  with  the  proviso 
that  instead  of  G(t»>)  we  read  G^(W), 

Next,  let  us  consider  the  case  where  the  input  to  the  system  of 
Fig.  1 is 

u(t)  s y(t)  [l4-m  g(t)J  4-  z(t),  (9.4) 

and  let  us  compute  the  cross-correlation  of  the  outputs  before  the 
luultiplier , which  we  denote  by  Pjnn4b^*^^*  notation  of 


cr 


Pmn*  = [ ih  f ) 

I ^ 


(9.5) 
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Using  the  methods  of  Sec.  5j  it  may  be  shown  that 


f r\  C.\ 
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x'lit!  x'etiCioj.'  SiiOllXd  Ou36rVS  "tlriGu  'tiiS  difl6r6nC6 


[p>intk  ~ ] 


(9.7) 


is  giverij  Up  io  a constantj  by  the  very  last  term  of  Eq.  (9»d). 
under  the  assumption  of  Eq.  (6.5).  Thus  the  average  mean  level 
out  of  the  system  due  to  the  presence  of  the  modulating  function 


is  increased  by  this  term  when  background  noise  is  present. 


However,  the  fluctuation  term  in  the  denominator  of  the  detection 


criterion  is  greatly  increased  by  the  background  noise.  Conse- 


quently, as  one  naturally  expects,  the  detectability  (as  given  by 
the  detection  criterjon)  decreases  as  the  background  noise  In- 
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10.  CONCLUSIONS 


In  the  foregoing,  we  have  developed  formulae  lifhlch  enable  us 
to  study  the  behavior,  and  the  optimum  design,  of  the  detection 
system  of  Fig.  1.  V^e  have  made  three  calculations  which  may  be 
of  general  interest.  We  have  (1)  introduced,  and  uerlved  various 
limiting  forms  for  the  deleotlon  crrrorioii,  (2)  computed  the  cross- 
correlation of  Che  outputs  entering  the  multiplier  for  two  kinds 
of  input,  stationary  Gaussian  noise  and  amplitude  modulated  noise; 
and  (3)  computed  the  autocori-elatlon  of  the  output  of  the  multiplier 
when  the  input  is  noise.  From  the  correlation  function,  the  corres- 
ponding power  spectrum  can  be  obtained  by  means  of  the  Wiener- 
Khintchine  relations - 

The  mathematical  techniques  used  here  may  be  of  use  in  many 
other  contexts  than  the  one  we  have  explicitly  considered.  In  the 
calculation  of  (2)  and  (3))  a basic  role  was  played  by  Eq.  (A?)  oi 
the  Appendix,  which  gives  an  explicit  expansion  of  the  higher  order 
statistics  of  a Gaussian  random  process  in  terras  of  its  second 
order  statistics.  By  using  this  expansion,  we  were  able  to  avoid 
using  the  Fourier  representation  of  a c;aiissian  random  process  used 
by  many  authors  ’ * It  appears  to  us  that,  as  long  as  only 

linear  and  quadratic  devices  are  considered,  the  use  of  the  Fourier 
representation  renders  many  computations  unnecessarily  cumbersome, 
and  may  not  always  readily  yield  the  correct  result  in  complicated 
situations  in  which  delta  functions  ai^e  involved  in  the  power 
spectra.  It  has  often  been  observed  that  the  correlation  function 
is  generally  better  behaved  than  the  power  spectrum  and  consequently 
the  mathematical  analysis  may  sometimes  be  simpler  if  one  first 
computes  the  correlation  function,  instead  of  the  power  spectra. 

The  computation  of  the  correlation  function  can  in  turn  be  facili- 
tated by  use  of  formula  (A7)» 


3 ■*  ...  .'u 


APPENDIX:  THE  MATHEMATICS  OF  NOISE 

In  this  appendix,  vie  siaminariy.e  the  main  notions  in  regard  to 
random  time  functions,  and  define  ’.vhat  is  mathematically  meant  by 
noise  ^ 

Let  y(t)  he  a fiinction  of  time  t defined,  for  the  sake  of 
generality,  for  - 4^  t < o»  - If  y(t)  contains,  in  addition  to  t, 

random  parameters,  so  that  at  each  point  t,  y(t)  does  not  have  a 
definite  value  hut  only  an  ensemble  of  possible  values  v/hich  it 
assumes  in  accordance  with  some  given  probability  distribution, 
then  y(t)  is  said  to  be  a random  time  function.  Then  the  ensemble 
average,  denoted  by  < > , is  an  average  taken  with  respect  to 

these  random  parameters.  We  let 

<Ai) 


r i-k, 

• u V V ^ < n •» 


(A2) 


denote  the  one^  and  two-point  ensemble  averages  of  y(t);  (t)  is 

called  the  mean-value  fxinction,  and  P„  (tj^,  t2)  is  called  the 
covariance  function  of  y(t). 

A random  time  function  for  which  the  ensemble  average 
^ y(t)  yCt  -P't  )y  has  a value  that  is  independent  of  t,  and  de- 
pends only  on  Tj  , is  called  stationary. 

The  autocorrelation  Rye'll  ) of  a random  time  function  is  de- 
fined (for  — v7-  ) by 


Tf  j-T  S ^ 

~ r Pu  ( t -t  v-s') "j 
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I.  (i 

where  the  subscript  TAv  is  used  to  denote  a time  average,  defined 
as  in  Eq.  A3« 
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Clearly,  for  a stationary  function 


Ry(-^  ) = TyCt,  t -t-r  ) 


(A4) 


The  power  spectrum  density  function  G (cd)  is  related  to 

w 

R (U)  by  means  of  a Fourier  transformation  (//iener-Khintchine 

y . :< 

theorem,"' . 
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By  noise  is  gtueraliy  meant  a stationary  ra.-dom  time  function, 
y(p),  whose  distrib’UPion  is  Gaussian,  The  important  mathematical 
property  of  such  a random  function  is  that  its  complete  probability'’ 
distribution,  and  all  possible  moments,  are  fully  determined  once 
we  know  the  mean  value  function,  /*=-.^(t),  and  the  autocorrelation 

V 

function,  R ,(”w  ).  The  mean  value  function,  jJL  (t),  is  often  assumed 

y 

to  be  identically  zero.  For  such  noise,  all  odd  order  moments 
vanish,  and  the  even  order  moments  may  be  expressed  in  terms  of 
the  second  order  moments  by  means  of  the  following  useful  formula. 

Let  n be  an  even  integer,  and  let  t^  ^ ,t„  be  points  of  time, 

some  of  which  may  coincide.  Then 

<'iW-  w(V/>  = 21 


where  the  sum  is  taken  over  all  possible  ways  of  dividing  the  n points 
into  n/2  combinations  of  2 pairs.  The  number  of  terms  in  the  summa- 
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tion  is  equa?.  to  l*3**’(n-3)  (n-1).  Thus,  for  nc  4 there  are 
three  terms  in  the  sum,  and  for  n - 8 there  are  105  terms. 

It  is  interesting  to  observe  that  Eq.  (A7)  characterizes 
noise.  A zero  mean  stationary  random  time  function  is  Gaussian 
if,  and  only  if,  all  its  odd  moments  are  zero,  and  its  even 
moments  satisfy  Eq.  (A7)^« 

By  random  noise  is  generally  meant  noise  with  a power  spectr'om 
that  is  constant  up  to  quite  large  frequencies.  The  covariance 
function  of  random  noise  is  thus,  more  or  less,  a J function. 
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